Calculation of the parity and time reversal violating interaction in 225 RaO 
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The 10-electron generalized relativistic effective core potential and the corresponding correlation 
spin-orbital basis sets are generated for the Ra atom and the relativistic coupled cluster calculations 
for the RaO molecule are performed. The main goal of the study is to evaluate the T,P-odd 
parameter X characterized by the molecular electronic structure and corresponding to a "volume 
effect" in the interaction of the 225 Ra nucleus Schiff moment with electronic shells of RaO. Our final 
result for X( 225 RaO) is —7532 which is surprisingly close to that in 205 T1F but has a different sign. 
The obtained results are discussed and the quality of the calculations is analyzed. The value is of 
interest for a proposed experiment on RaO [PRA 77, 024501 (2008)] due to a very large expected 
Schiff moment of the 225 Ra nucleus. 



INTRODUCTION 



structure of the molecule: 



Study of fundamental interactions, which break both 
the time-reversal invariance (T) and spatial parity inver- 
sion symmetry (P), or T,P-odd interactions, is a way to 
study so-called "new physics" [l|, 0] beyond the Standard 
model of electroweak and strong interactions. Despite 
well known drawbacks and unresolved problems of the 
Standard model (radiative corrections to the Higgs mass 
are quadratically divergent; rather artificial Higgs mech- 
anism of symmetry breaking is not yet proved experi- 
mentally; the origins of CP-violation, where "C" is the 
charge conjugation symmetry, are not well understood; 
the CP-violation according to CPT-theorem is equivalent 
to T-violation problem, etc.) there are very few exper- 
imental data available which are in contradiction with 
this theory. 

In turn, some popular extensions of the Standard 
model, which allow one to overcome its disadvantages, 
are not confirmed experimentally. Since considerable en- 
hancement of the T,P- violation effects is expected in po- 
lar diatomic molecules comprising a heavy atom, such 
systems have been a research subject for a couple of 
decades. It is also believed that the Standard model 
gives the baryon-to-photon ratio 10 orders of magnitude 
smaller than the observed value. This is a strong argu- 
ment in favour of additional source of CP-violation which 
may be detected using atomic and molecular experiments 
searching for T,P-violating interactions. Thus, finding a 
molecule that exhibits the strongest T,P-odd effects may 
be crucial for observation of these effects experimentally. 

Following 3] , the effective interaction with the Ra nu- 
cleus Schiff moment in RaO can be written in the form 
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where (Tn is the Ra nuclear spin operator, A is the unit 
vector along the internuclear axis z from Ra to O, S is the 
Shiff moment of Ra, X is determined by the electronic 



x,y,z— 



where (r) is an electronic density calculated from the 
four-component wave function ip. The amplitude of the 
T,P-odd spin-axis interaction, 6SX, in RaO was esti- 
mated by Flambaum [4] to be 500 times larger than that 
in TIF [5J|. Such an enhancement in RaO js based firstly 
on the 225 Ra nucleus Schiff moment la, M| being about 
200 times greater than that of 205 T1 [8l4lOl| and secondly 
on a semiempirical estimate X( 225 RaO) = 2.2X( 205 T1F). 

In the paper we have performed ab initio calculations of 
the T,P-odd parameter X( 225 RaO), excitation energies 
of the Ra atom and some spectroscopic properties of the 
225 RaO molecule. 



METHODS 

It follows from eq. [2] that the volume effect is local- 
ized on the nucleus and determined by polarization of 
the valence electrons. Earlier it was demonstrated by our 
group [ll[ that calculation of such core properties can be 
performed efficiently in two steps. Firstly, electron cor- 
relation for valence and outer-core electrons is accounted 
for, whereas the inner-core electrons are excluded from 
this calculation using the generalized relativistic effective 
core potential (GRECP) method 12|,|l3(, which yields a 
very accurate valence region wave function by the most 
economical way. Secondly, since the heavy atom's inner 
core wave function cannot be accurately obtained within 
the GRECP method, it has to be recovered using a non- 
variational restoration procedur e II ill. The two-step ap- 
proach has recently been used in |l4l4l8j for calculation of 
other core properties, such as hyperfine constants, elec- 
tron electric dipole moment enhancement factor, etc., in 
molecules and atoms. 
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The fully-relativistic Fock-space coupled cluster code 
with single and double cluster amplitudes (FS-RCCSD) 

IFA 2pj was applied to take account of both the electron 
correlation and relativistic effects. Triple cluster ampli- 
tudes and basis set enlargement corrections to X were 
obtained using scalar-relativistic approach within CFOUR 

2l| and MRCC 22|, |23| codes via interface to our new 
non- variational one-center restoration codes developed in 

171. 



GRECP generation and atomic calculations 

The GRECP for Ra is generated to simulate interac- 
tions of ten explicitly treated electrons, the valence and 
outer-core (6s6p7s) ones, with the inner-core Is — 5d elec- 
trons, which are excluded from calculations. 

Basis sets for Ra were constructed using the general- 
ized correlated scheme [2^ : as a result, five correlation 
spin-orbital basis sets were generated, which are desig- 
nated here as A, B, C, D and E. Each one is constructed 
of optimized primitive 20s, 20p, lOd, 8/ and 5g Gaussian- 
type functions contracted to [4s4p3dl/lg], [5s6p4d3flg], 
[6s7p5dAf2g], [6s8p4d2flg] and [6s8p5d5f3g] sets, re- 
spectively. These basis sets were obtained in a series of 
atomic two-component GRECP calculations, with cor- 
relations included by the FS-RCCSD method 0, and 
optimized to reproduce the 7s 2 — > 7s7p and 7s 2 — > 7s6d 
transition energies of the atom as accurately as possible. 

Our FS-RCCSD results are compared to experimental 
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26j in the table Q] A reasonable level of accuracy 
has been achieved with C and E basis sets. The use of D 
basis set leads to satisfactory results as far as 7s 2 — > 7s7p 
transition energies are concerned. However, a lack of / 
and g functions apparently causes considerable errors in 
7s 2 — > 7s6d transition energies. 



TABLE I. Transition energies (TE) from GRECP/FS- 
RCCSD calculations of the lowest-lying states of the Ra atom 
for 10 correlated electrons. All values are in cm - . 



7s 
7s 
7s 
7s 
7s 
7s 
7s 
7s 



/2 7Pl/2 
/2^P\/2 
/27P3/2 
/2^P3/2 

/2®d, 5 / 2 
/2®d 5 / 2 



(3=0) 
(3=1) 
(3=1) 
(3=2) 
(3=1) 
(3=2) 
(3=2) 
(3=3) 



12435 
13386 
21546 
16262 
15579 
15969 
19864 
16901 



12920 
13840 
21258 
16660 
14112 
14434 
17943 
15191 



13072 
13995 
21045 
16756 
13892 
14199 
17553 
14925 



13164 
14136 
21714 
16996 
14964 
15318 
18946 
16207 



13029 
13953 
21000 
16714 
13753 
14049 
17340 
14761 



13078 
13999 
20715 
16689 
13716 
13994 
17081 
14707 



Molecular calculations 



To perform relativistic Fock-Space coupled cluster cal- 
culation of RaO, the basis set D was used instead of E, 
thus, excluding several d, f and g basis functions to re- 
duce the computational efforts. Note that the major con- 
tribution to the X value is determined by the s— p orbital 
mixing, therefore, such simplification is justified here. 
The (10s5p2dlf)/[As3p2dlf] basis set from the molcas 
4.1 library was used for oxygen. 

A one-component self-consistent-field (SCF) calcula- 
tion of the (la . . . 5cr) 10 (l7r27r) 8 ground state of RaO 
is performed first. It is followed by two-component 
RCC calculations taking account of single and double 
(RCCSD) cluster amplitudes. To estimate the conver- 
gence of the results with respect to the basis set enlarge- 
ment we have performed two scalar-relativistic calcula- 
tions using the coupled cluster method with single and 
double cluster amplitudes in the basis set D and in en- 
larged basis set (uncontracted basis set D) using CFOUR 



code [21 1 . 



Correction for the contribution of the triple cluster 
amplitudes was also estimated. For this we have per- 
formed (i) scalar-relativistic calculation using the cou- 
pled clusters with single and double cluster amplitudes 
and (ii) scalar-relativistic calculation using the coupled 
clusters with single, double and triple cluster amplitudes 
(CCSDT). The final X value is obtained as 



A (FINAL) = A(RCCSD) + A(CCSD uncon.) 
+ A(CCSDT) - 2A(CCSD). 



(3) 



The GRECP/RCCSD method was also used to cal- 
culate the ground state potential curve of the RaO 
molecule, which was then used to obtain spectroscopic 
properties of RaO. The results are listed in table Hill All 
molecular calculations of X were carried out for the equi- 
librium internuclear distance, R e — 3.852, given also in 
table [OH 



Leading 






Basis set 


Exper. 




conf. 


Term 


A 


BCD 


E [25, 26] 


TABLE II. 


7s 1/2 -> 










state R e = 



Parameter X calculated for the 225 RaO ground 
3.852 a.u. Only the GRECP/RCCSD result takes 
into account the spin-orbit interaction. All values are in a.u. 



Method 



X 



Hartree-Fock -9609 
RCCSD -7696 
CCSD -7648 
CCSDT -7209 
CCSD uncontracted -7923 
FINAL -7532 
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TABLE III. Calculated properties of RaO: equilibrium inter- 
nuclear distance (R e ), harmonic frequency (w e ) and vibra- 
tional anharmonicity (w e x e ). 



Re, A 


3.852 


w e , cm -1 


598 


w e x e , cm -1 


-1.09 



than that expected from the speculative RaO experi- 
ment. However, considering the impressive progress in 
molecular spectroscopy during last decades one can ex- 
pect an order of magnitude or even better sensitivity of 
the RaO experiment than that attained on TIF. 

Analysis of contributions to X 



RESULTS AND DISCUSSIONS 

Although triple and higher cluster amplitudes in the 
valence region are important for chemical and spectro- 
scopic properties in general, it can be concluded from 
table HI] that GRECP/RCCSD calculation includes all 
major correlation contributions to X, since triple clus- 
ter amplitudes contribute only about 6%. Even smaller 
correction was found for basis set enlargement. We do 
not expect that further enlargement of the basis set and 
accounting for quadrupole amplitudes will change the re- 
sult by more then 1%. The influence of the inner core — 
valence electron correlation on X value was estimated in 
27| to be no more than 2% for TIF. Taking into account 
essentially different electronic structures of TIF and RaO 
we believe that our final result for X is reliably valid 
within 10% of accuracy. 

One of the goals of this work is to determine whether 
or not the RaO molecule is a feasible candidate for ex- 
perimental research of T,P-odd effects. The corrected 
value X( 225 RaO) « -AT( 205 T1F), according to our cal- 
culations, is not as large as could be expected from the 
usual Z-scaling (see [28| and references). The estimation 
based on Z-scaling determines only a "systematic" vari- 
ation with increasing Z, and does not take into account 
some particular changes from one element to the next 
in the Periodic table and features of chemical bonding 
with different elements (F in the case of TIF and O in 
the case of RaO). A very large estimated Schiff moment 
of 225 Ra @,0 will strengthen the T,P-odd interactions 
in RaO by two orders of magnitude. The best upper 
bound on the nuclear T,P-odd interactions can be ob- 
tained using the limit on the 199 Hg EDM, |d( 199 Hg)| < 

3.1 x 10~ 29 e-cm[2!|, and the atomic calculation that links 
the atomic EDM with the Schiff moment of 199 Hg nucleus 
[13, d( 199 Hg) = -2.8xl0~ 17 (S( 199 Hg/efm 3 )e-cm. Com- 
bining these values, the current upper limit on |S( 199 Hg)| 
is 1.1 x lCP 12 /efm 3 . Using the corresponding experi- 
mental datum for TIF [3(| and the most accurate calcu- 
lated value for X Q one can obtain the upper limit on 
|S( 205 T1)| as 1.7 x lCT 10 /efm 3 . Supposing the same ac- 
curacy for a "speculative" RaO experiment as for TIF 
(performed twenty two years ago), the same limit on 
|S( 225 Ra)| will be obtained. However, taking into ac- 
count the relation S( 225 Ra)/S( 199 Hg) « 200 [E-E^ and 
the upper bound on |S( 199 Hg)| we obtain |S( 225 Ra)| < 

2.2 x 10~ 10 /efm 3 which is only slightly (1.3 times) worse 



In the present paper an analysis of contributions to X 
based on the one-center representation of the one-particle 
density matrix has been performed assuming the scalar- 
rclativistic approximation used in the paper to construct 
initial (reference) Hartree-Fock wavefunction and molec- 
ular orbitals. The analysis is similar to that performed by 
our group for Eu 2+ in Ref. 17 1. It should be noted that 



such one-center density matrix analysis is free from those 
rotations in the space of molecular orbitals which do not 
change the wavefunction of a considered state (e.g., rota- 
tions between the doubly occupied orbitals) and can also 
be easily performed in the case of a correlated wavefunc- 
tion, i.e. beyond the one-configuration approximation. A 
mean value of the one-electron operator X can be eval- 
uated as follows: 



pqXqp 



(4) 



where X pq are matrix elements of X, D pq are matrix 
elements of the one-particle density matrix, and indices 
p, q enumerate all the atomic basis functions centred on 
Ra and O assuming the conventional MO LCAO approx- 
imation used in the RaO calculation. For a qualitative 
analysis of contributions to X it is reasonable to rewrite 
Eq. (j4j in terms of only atomic orbitals of Ra, both occu- 
pied and virtual, which can be calculated using Hartree- 
Fock method for a single Ra atom since such an analysis 
can be performed excluding one more ambiguity caused 
by overlapping of the atomic orbitals on different centers. 

Represent the molecular density matrix D pq in a "suf- 
ficiently complete" basis set of these atomic orbitals with 
indices i, j, running over all the occupied and virtual or- 
bitals of Ra , i.e., decompose all original basis functions 
p, q (centered either on Ra or O) into one-center func- 
tions Then the mean value for X can be rewritten 
in the form 



(X) 



z2 /; "- v " 



(5) 



where Xij arc matrix elements of X in the enlarged basis 
set of Ra orbitals, described above. (X) can be slightly 
different from (X) due to incompleteness of the used en- 
larged Ra basis set. In practice however, this difference 
can be made negligible. 

Applying this scheme to RaO we obtain the following 
contributions: the main one (more than 50%) is due to 
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a polarization of outer-core 6p z orbital of Ra (z-axis co- 
incides with the axis of the RaO molecule) into 7s, 8s.... 
The rest are due to mixing of 7s, 8s, 9s . . . orbitals of Ra 
with different p z orbitals of Ra which partially can con- 
sist of orbitals of oxygen (in terms of the basis set p, q 
used in RaO calculations). 

It should be noted that this kind of analysis can not 
be performed directly for individual contributions from 
canonical molecular orbitals (see table IIV[) due to the 
above mentioned large unitary rotations between atomic 
orbitals of radium (leading to mutually compensated un- 
physical contributions) and their overlapping with or- 
bitals of oxygen. These circumstances result in the pres- 
ence of several leading contributions. This is opposite to 
the case of TIF (see Ref. @) where, due to the relatively 
simple electronic structure, such orbital analysis was suc- 
cessfully applied with only one leading contribution. 



and X(RaO) discussed earlier we can write 

X(T1F) = 2 < 6s\V\6p >< 6s\W\6p > /AEdf, (6) 
AT(RaO) = 2 < 6p\V\7s >< 6p\W\7s > /A£ Ra0 , (7) 

where AE TW = E 6s - E 6p , AE R&0 = E 6p - E 7s , V is 
the polarization operator, W is an operator for the X 
property localized on the Ra and Tl nuclei (see Eq. ((5J)). 
Note that AStif and AE^o have the same signs. Let 
us choose the phase of orbitals so that the first ex- 
ternum would have a positive sign for all the functions 
of Ra and Tl. Then the matrix elements < 6s|14 y |6p > 
and < 6p|W| 7s >, localized on the nucleus, will have 
the same signs. However, the terms < 6s|y|6p > 
and < 6p|y|7s >, localized in the valence region, will 
have different signs, which explains the opposite signs of 
X(T1F) and X(RaO). 



TABLE IV. Parameter X calculated for the 225 RaO ground 
state R e = 3.852 a.u. Individual shell contributions are cal- 
culated from the spin-averaged GRECP/SCF orbitals. All 
values are in a.u. 



Shell : main contribution 

6s 2 (Ra) 
2s 2 (0) 
6pl(Ra) 
2pl(0) 
6p 2 6p 2 (i?a) 



la' 
2a 2 
3a 2 
4a 2 
5a 2 
Itt 4 

2tt 4 : 2pl2p 2 (0) 
Total SCF(spin-averaged) 



X 



-50 
3862 
9936 
-14426 
-9036 
50 
54 
-9609 



Analysis of X sign 

The different sign of the X in TIF and RaO molecules 
can be explained in a simple qualitative analysis. The va- 
lence electronic configurations of Ra and Tl are 6s 2 6p 6 7s 2 
and 6s 2 6p 1 , respectively. TIF has ionic bonding, the 6p 
valence orbital of thallium becomes mainly unoccupied 
due to interaction with fluorine. Thus, qualitatively, the 
electronic configurations of the TIF molecule can be writ- 
ten as Tl+(6s 2 )F-(2s 2 2p 6 ). Electronic structure of RaO 
is more complicated. Calculation shows that its config- 
uration is rather close to Ra + (6s 2 6p 6 )(7r 4 )0 - (2s 2 2p 2 ), 
where n is three quarters 2p x i y \[0) and one quarter 
6dxz(yz) (Ra) ■ According to the Mulliken population anal- 
ysis, an effective configuration of Ra in RaO is 6s 2 6p 6 6<i 1 . 
Thus, due to the interaction with oxygen, the 7s valence 
orbital of Ra becomes mainly unoccupied in RaO. Follow- 
ing the density matrix analysis the leading contribution 
to X in TIF and RaO is due to the polarization of 6s and 
6p orbitals, correspondingly (see above). 

Taking account of the leading contributions to X(T1F) 
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